Engineering Superposition States and Tailored Probes for Nano-resonators Via 

Open-Loop Control 
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We show that a nano-resonator can be prepared in mesoscopic superposition states merely by 
monitoring a qubit coupled to the square of the resonator's position. This works for thermal initial 
states, and does not require a third-order nonlinearity. The required coupling can be generated 
using a simple open-loop control protocol, obtained with optimal control theory. We simulate the 
complete preparation process, including environmental noise. Our results indicate the power of 
open-loop control for state engineering and measurement in quantum nano-systems. 
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The quantum behavior of nano-mechanical resonators 
is an area of increasing activity [H [21 [3l IH |5]. These 
devices offer the exciting prospect of observing quantum 
behavior in a macroscopic mechanical system [SJ |7] . They 
can also be integrated with superconducting circuits, and 
have potential applications such as sensing [B^ and in- 
formation processing [9]. In addition, superconducting 
"stripline" resonators are now realizing quantum effects 
previously the exclusive domain of quantum optics. Re- 
cently two major advances have been made in this area: 
observation of the discrete energy levels of a stripline res- 
onator [To] , and a controllable single-photon source [TT] . 

The task of preparing a nano-mechanical resonator in 
a nonclassical state has yet to be achieved experimen- 
tally. Here we consider preparing a superposition of two 
spatially separated wave-packets, a primary goal in this 
field. For a resonator to evolve to such a state itself re- 
quires that it possess a third-order nonlinearity {{a^a)'^ or 
equivalently an a;^ potential) [12J, and the nonlinearities 
that exist in nanoresonators, at least at present, are too 
weak for this purpose [T^]. A technique presented in [14 
allows a resonator to be prepared in an arbitrary state via 
an auxiliary qubit, but this is best suited to states with 
only a few phonons. Here we show that superposition 
states can be prepared using an entirely different process. 
This requires only a continuous measurement of a qubit 
that is sensitive to the square of the resonator's position. 
We show that this preparation process can implemented 
by applying open-loop control to the qubit |3I]. This 
approach was inspired by [TS], where it was shown that 
control sequences applied to an auxiliary system could be 
used to generate nonlinearities in a target system. While 
the work in |15j provides a proof of principle, we sug- 
gest that a better approach is to find control protocols 
using "optimal control theory" (the name given to a spe- 
cific optimization technique for open-loop control proto- 
cols) [T^ . This has been exploited to great effect in fem- 
tosecond control of chemical interactions [17j , and more 
recently for constructing quantum gates and entangling 
operations in mesoscopic systems fTE\ [TOl [50] . The pri- 



mary question is the time-scale on which the control must 
operate, as this determines the feasibility of the proce- 
dure. While we explicitly consider mechanical resonators 
in what follows, the techniques apply, virtually without 
alteration, to superconducting stripline resonators. 

We now turn to our first result, that mesoscopic super- 
position states can be prepared purely by monitoring the 
square of a resonator's position. Under such a measure- 
ment, in addition to the evolution due to the resonators 
Hamiltonian, the density matrix, p, evolves as 21j 

dp = -k[x^, [x^, p]]dt + V2k{{x^, p} - 2{x^)p)dW (1) 

where {x^,p} — x'^p + px^ denotes the anticommutator, 
X = a-|-a^ is the dimensionless position operator (a is the 
resonator annihilation operator), and k is the strength 
of the measurement (characterizing the rate which the 
measurement extracts information). 

To provide a physical realization of a measurement of 
x^ we employ a Cooper-pair box (CPB) as a probe sys- 
tem: we couple the CPB to the resonator via the interac- 
tion pozx"^ (we will consider how to engineer this coupling 
below) , and perform a continuous measurement of Ux on 
the CPB with strength n. This measurement is described 
by Eq.(fl]), with x"^ replaced by a^, and can be realized 
by coupling the CPB to a superconducting stripline res- 
onator, a method that has already been demonstrated 
with excellent resolution [TT]. The analysis in [23] in- 
dicates that K ~ 10^ s~^ is realistic, given an average 
value of 10^ photons in the stripline. Here we choose 
K = 8p — 8w', where uj' is the effective frequency of the 
resonator as seen by the CPB [35j. In all our simulations 
we choose uj' = 27r/' = 27r//100, where / = 100 MHz is 
the resonator frequency. 

For reasons that will be explained later, we implement 
the indirect measurement of x^ by alternating between 
the interaction and the measurement of a^, with each 
switched on for a duration of 1/(160/). We simulate 
this process with the resonator at both zero tempera- 
ture and T = 22.7 mK, with the resonator quality fac- 
tor Q = 10^ [35]. At this temperature resonator's mean 
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FIG. 1: (Color online) Snapshots of the Wigner function for 
a harmonic oscillator with effective frequency /' = 1/r = 
IMHz, subjected to a continuous measurement of the square 
of position. Luminosity denotes the absolute value of the 
Wigner function (online: blue is positive and green negative). 
The X-axis is dimensionless position, x — a + a\ and the y- 
axisisp = — j(a — a^). Both axes cover the interval [—7.1,7.1]. 
Top row, temperature T = 0: a) t = (vacuum state); b) 
t = 6.25r; c) t = 13.75r. Bottom row, T = 22.7 mK: d) i = 
(thermal state); e) t = 17.625r; f) i = 18. 5r. 



phonon number is (n) = 4.2 [29 . These numerically 
intensive simulations were performed using a regenera- 
tive Monte Carlo algorithm recently developed by one of 
us, running on a 128-node parallel computer 34J. Snap- 
shots of the Wigner function for the resonator in these 
two cases are shown in Fig. [T| At T = a coherent su- 
perposition of two localized wave-packets is formed, and 
these slowly increase in separation. At T ~ 23 mK this 
superposition is still created, although it looses its co- 
herence when the separation is sufficient for the ther- 
mal decoherence to overpower the purification due to the 
measurement. The coherence can be preserved for larger 
separations by increasing either Q, ot fi and k. 

To understand the creation of superposition states, 
first note that a measurement of x'^ does not distinguish 
between between positive and negative x. Thus oscillat- 
ing states whose probability densities are symmetric in 
position are stable under this measurement. The super- 
position states that are created have this form. Second, 
whenever a continuous measurement of x^ produces a 
non-zero result i^ for any appreciable time on an initial 
vacuum state, which it must do at some point because the 
stream of measurement results is random, this produces 
a state with peaks close to x = ±x. This is because i^ 
becomes the most likely square displacement of the res- 
onator. The resulting symmetric superposition tends to 
persist because it is stable under the measurement. This 
process continues, in which new measurement results ran- 
domly change the twin peaks of the position density. The 
separation of the peaks thus undergoes a random walk on 



the real line, with the addition of a slow drift upwards. 
This drift is due to the fact that the measured observable 
does not commute with the Hamiltonian, and so feeds 
energy into the system. This can create superposition 
states with arbitrarily large separation. 

To generate superposition states using the above 
method, we must realize the coupling ^azX^ between the 
resonator and the CPB. One can readily obtain the in- 
teraction a^x merely by placing a charge on a resonator 
situated adjacent to the CPB. To obtain an x^ interac- 
tion, we note that it was shown in |15| that, given the in- 
teraction <TzA between an auxiliary qubit and a "target" 
system, an effective Hamiltonians for the target of the 
form j4" could be generated by manipulating the qubit. 
The method in [TS] can be adapted for the present prob- 
lem, but it is not especially practical. 

To provide a realistic control protocol for use in a solid- 
state circuit, given the constraints on switching rates or 
other resources, we apply optimal control theory [16 . In 
our case we want to generate evolution corresponding to 
the effective Hamiltonian Hcs — hw' a) a + hfiazX^, given 
the physical Hamiltonian 

Hphys = huja* a + hA{t)<TzX + AEax- (2) 

Here AE is the Josephson energy of the CPB, and A 
is the strength of the resonator-CPB interaction. The 
CPB is operated at the degeneracy point where the qubit 
is largely protected from the low- frequency noise. We 
first modulate the interaction strength so that A(t) = 
2Acos(t^i), with 1/ — AE/h — {lo — oj'). Moving into the 
interaction picture with respect to the CPB, and drop- 
ping the high frequency terms, results in the Hamiltonian 



. + hXazX + Hcti-\{t), 



(3) 



where we have now added a time-dependent control 
Hamiltonian. We take the control Hamiltonian to be of 
the form -ffctri = !T'[cx(t)ax + Cy{t)ay]. This Hamiltonian 
can be implemented by applying a time-dependent mag- 
netic field to the SQUID loop of the qubit [22] . We expect 
that for the control protocol to work well, the physical 
interaction strength should be significantly larger than /i. 
We therefore chose (arbitrarily) A = 200/i — 200a;'. With 
uj' — 27rMHz, this means A = 47rx 10* s~^, which is about 
a factor of 10 higher than the estimate for realistic values 
given in [53] for nano-mechanical resonators, but realis- 
tic for coupling with superconducting resonators |10l lllj . 
We then use optimal control theory to search for a con- 
trol Hamiltonian to generate iJeff with fi ^ lu' . Specif- 
ically, we search for functions Ca,(i) and Cy(t) to gen- 
erate the evolution exp(— i_ffeffAt), for a time interval 
At = 1/(320/): this value is chosen so that At < 1//, 
to realize a quasi-continuous measurement. Note that the 
larger A and k are, the larger is the measurement rate for 
x^, important to beat thermal noise. 

The control Hamiltonian we have obtained is shown in 
Fig. [2] For the majority of the protocol, the rate at which 
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FIG. 2: A control protocol for generating the effective inter- 
action fia-^x^ between a nano-resonator and a qubit, where 
the physical interaction is 200 fia^x, and the effective fre- 
quency of the resonator is lu' = fi. The control Hamiltonian 
is H{t) — h[cx{t)ax + Cy{t)ay], with Cx shown in (a), and Cy 
shown in (b). The insets give an expanded view of the first 
tenth of the protocol. 



Cx and Cy must be varied is less than 100/i per Ai/10. 
(The exception to this is in the first 40th of the protocol, 
where the Hamiltonian is required to change about ten 
times faster.) For /x = 27r x IMhz, this gives a rate of 
change of IGHz per 3ns, which is readily achievable with 
current CPB circuits jUJ US]. This is probably about 
the largest interaction strength that one can reasonably 
engineer with this protocol; taking fj, = 2tt x lOMHz re- 
quires a rate of lOGHz per 0.3ns. This can be done with 
fast electronics, but would be restricted by the induced 
voltage in the SQUID loop [3S] 

It turns out that the protocol presented above does not 
generate Hcs with quite enough accuracy. To fix this we 
modify the protocol slightly to improve its symmetry: we 
apply it once, flip the sign of the resonator-CPB interac- 
tion, and then apply the protocol again with the signs of 
Cx and Cy flipped. This symmetrizes the evolution gener- 
ated by Hcs under a flip of the CPB charge states. The 
resulting protocol takes a time of 1/(160/), plus the time 
to flip the sign of the interaction. 

To generate the effective x'^ measurement, we rapidly 
switch between the control protocol and a continuous 
measurement of the x-component of the CPB, at time in- 
tervals of 1/(160/), just as we did previously. The reason 
we alternate the measurement with the control protocol 
is that the measurement interferes with this protocol if 
they are performed together. 

We flrst simulate the entire open- loop control process, 
with the resonator at zero temperature, and display a 
snapshot of the resulting evolution of the nanoresonator 
in Fig. p^ (a). The procedure produces a superposition 




FIG. 3: (Color online) The Wigner function for a harmonic 
oscillator measured via a Cooper-pair box (CPB). Plots (a)- 
(c): Mesoscopic superposition states generated by the mea- 
surement for various CPB decoherence rates, 7, and T = 0. 
The phase-space scale is the same as Fig. [I] (a) 7 = 0; 
(b) 7 = 0.1/'; 7 = 0.5/'. Plots (d)-(f) are the respective 
side-views of plots (a)-(c). Plots (g)-(i) are a sequence where 
the resonator is initially prepared in the ground state, with 
7 = 0.1/' and the bath at T = 30mK. 



state as expected, although this state is no longer sym- 
metric, (the heights of the two wave-packets are differ- 
ent). This asymmetry is not especially important, so long 
as both packets are appreciable. In Fig. |3jd) we show a 
side view of the Wigner function, giving the heights of 
the interference fringes relative to the wave-packets. 

To complete our analysis we consider the combined ef- 
fects of environmental dephasing on the CPB, and non- 
zero temperature. The decoherence rate due to this de- 
phasing is much faster than the thermal damping rate of 
the resonator. To date, the best realizations of CPB's 
have a decoherence rate of about lO^s^^ at the degener- 
acy point |26l 127] . It is therefore important to examine 
the effects of this decoherence on our control protocol. 
To this end we simulate the entire measurement/control 
process, including dephasing of the CPB qubit in all 
directions at rate 7. This requires adding the term 
— 7^Jcri, [ci, p]]dt, where i = x,y,z, to the evolution 
of the density matrix. We choose dephasing equally in 
all directions to represent a worst-case scenario. 

We flrst simulate the evolution when the resonator is 
at zero temperature. The results are presented in Fig. [3] 
(b/e) and (c/f), for 7 = lO^s'^ and 7 = 0.5 x lO^s"!, 
respectively, with /i = 27r x IMHz as above. We see 
that this rate of decoherence does effect the superpo- 
sition states generated by the measurement, but does 
not destroy the superposition completely; the interfer- 
ence fringes still exist, albeit signiflcantly reduced. We 
then simulate the evolution when the resonator is at 



30 mK, with a CPB decoherence rate 7 = 10^s~^. In 
this case we prepare the resonator in the ground state 
before we start the measurement process. A number of 
cooling schemes have been proposed that can potentially 
perform this preparation (see, e.g [3D]). Superposition 
states can still be created in this case, although they do 
not last long, as is shown by the sequence of snapshots 
in Fig. [3] plots (g)-(i). These snapshots are separated by 
a quarter of an oscillation period. 

To summarize, we have shown that a measurement of 
the square of a resonators position will generate meso- 
scopic superposition states, even when the resonator is at 
finite temperature. We have also presented a method to 
realize this measurement, using open- loop control. While 
our example is probably the simplest nonlinear interac- 
tion to create in this way, our results indicate that open- 
loop control protocols, derived using optimal control the- 
ory, may have considerable potential for engineering a 
range of non-linear interactions and measurements for 
quantum systems. While we have used a Cooper-pair 
box to couple to the resonator, there are other systems 
that can be used for this purpose. Polar molecules, for 
instance, can be coupled to nano-resonators in essentially 
the same way [28], and suffer far less decoherence. 

The main limitation on our open-loop control proto- 
col is the time-scale required for the control. While the 
protocol we have presented here is adequate for our pur- 
poses, an important open question is whether there exists 
a protocol for this task that allows slower control frequen- 
cies. More generally, one would like to place lower bounds 
on the control frequencies required for many more tasks 
relevant to quantum state-engineering and measurement, 
and this will be the subject of future work. 
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